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TITLE: Hierarchial Analysis of Variance 



AUTHOR: Richard A. Lamm 



ABSTRACT : 



This routine performs a complete analysis of variance for data from a 
repeated subsampling design with an arbitrary number of stratifications 
and arbitrary number of observations within the stratifications. The 
amount of data analysed is restricted only by the available drum memory. 



DISCLAIMER: 



"The authors of this program material, the POOL organization 
and Royal McBee believe this program to be correct; however, 
they bear no responsibility, financial or otherwise, for 
errors resulting from its use. This program is distributed 
only to individual and installation members of POOL. Further 
distribution of this manual and accompanying tapes for use 
by non-members is prohibited." 



FOR POOL MEMBERS ONLY 



Issued March 20, 1961 Printed in U. S. 



Program No.- Fh-l65 

1. Title: Hierarchial Analysis of Variance 
Author: Richard A. Lamm 

Date: December 7, 1959 

Installation: Lederle Laboratories, Pearl River, Pew York 

Classification:. F-U 

2. Purpose: 

To perform a complete analysis of variance for data from a repeated 
sub sampling (nested subsampling) design with an arbitrary number of 
stratifications and arbitrary numbers of observations within the 
stratifications. Transformation of the data to 1000 log x is control- 
led by the Transfer Control button. 

3. Restrictions: 

a. All observation data must be entered (coded) at a q of 30. Generally, 
data should be coded so that they are k decimal digits in length. 

b. No (coded) Sum of Squares for a level may exceed 230 _]_ # 

c. No (coded) mean may exceed 2-*-° _1. 

d. The amount of data is restricted only by the available drum memory. 

e. The number of lines in the Analysis of Variance is limited to 11. If 
more lines are required, it is necessary to change L0+O6I4.2 from XC6300 
to XCttss in order to provide 3L consecutive sectors for storage of the 
analysis, where L is the number of lines. 

f. There must be at least two observations within each cell of the lowest 
heirarchial level, since the Zx and Zn^-j^... are stored sequentially 

in pairs in the location of the data. This restriction is not necessary 
if there is a sufficient number of observations in the beginning cells 
to prevent the storage of partial sums on unused data. 
However, a 'divide check 1 will occur at Lo / 22 of FI1-I6I4 if n is 1. 
Press start compute to continue. The result will be correct. 

g. Memory Allocation. (Routines in parenthesis are not used with this 

program) 

0000-0263 (Program Input, 10. k) 

0300-0563 Data Input No. 3, 11.2 

0600-0963 (Decimal Memory Printout) 

1000-1058 Integer Printout 

1100-1231 (Data Output No. 1, 12.0) 

1232-1263 Modifications to Data Input No. 3 

1300-1531 Data Output No. 2, 12.1 

1532-1563 Storage of Data Input No. 3 Modification Results 

1600-1663 (Square Root, 15.1) 

1700-1863 Log k x, 18.0 

I9QO-I963 Alphanumeric, 19.0 

2000-2063 (Modifiable Hexadecimal Input) 

2100-2163 (Antilog Printout) 

2200-2263 Zero Suppression for 12.0 and 12.1 

2300-2363 (Exponential, 17.0) 

2l|00-2li63 (Antilog) 

2500-2663 Data Output "30" 

2700-5963 Available for Program and Data Storage 

6000-6263 Double Precision Sum of Squares or Products 

6300-6332 Analysis of Variance Results 

6333-6363 Intermediate Results-See Attached Code Sheet for Track 63 



I/! tie: ?:iere.rchlal Analysis of Variance program N". - :'4-lO:.: 

December ?, 1?5'9 

h. If it Is desirable to obtain the irjean square for each of the in- 
dividual level 1 cells (level Mean Square), the following char,,.;.. 
are required; 
Location Present Charge To 



001+7 Y012U Y0507 

022k A[ ] UC460 

0152 B0120 U0509 
(All addresses are understood to be Lo+) 

Method: 

The method used is the same as may be found in many books on 
statistical methods, e.g., R. L. Anderson and T. A. Bancroft. 
Statistical Theory in Research . McGraw-Hill Book Company, Inc., New 
York (1952;. See paragraphs 22.5 and 22.4. 
The sums of squares are obtained as follows: 

«. y yyfTy 1 - (ST 

L . Z_, A- ^— ^N T>,;k...< 



Level 1 



etc. 



s j 1 L 1 "^tjK.-S 

y yy(2if_y y (jjgf 

V-VV ' Tli i K -' s V"T ^.m ti K...s 



5. Coding Information: 

a. Storage. The program requires 7 tracks of storage for instructions 

and constants. In addition, it uses 63OO to 63OO+3L-I, 
where L is the number of lines in the Analysis of Variance 
and 6333, (651*2), 63W, (6346), 63^7, 63^9, 6350, (6355), 
6356-6O, 6362 and 6363. The addresses in parenthesis are 
used only with the alternate output. 

b. This program uses: 

Subroutine Location 

Data Input No. 3 with Lederle Modifications 0300 

Integer Printout 100 ° 

Data Output No. 2 1300 

Logfc x (If transformation is desired) 1700 

Alphanumeric 1 9°° 

Data Output "30" 25OO 

Double Precision Sum of Squares or Products 6000 

There is no calling sequence required for this program. Halt and 
transfer to Lo for computations. 

c. Input. 

A data tape should contain; 

1. Experiment designation or heading followed by a conditional stop. 

2. The coding of the data, e.g. 1000 x, followed by a conditional 
stop. There should be no carriage return preceeding this inform- 
ation. 
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3. Data as follows: 

a. 0+29(I^o d ) , nlll... , ... , niJ;k... , -0000000' 

where nijk... is the number of observations in the ith cell 

of the first level, in the jth ce n of the second level, in 

the kth C ell of the third level, etc. 

O+296363 *n.n. . . » . . . » n . Jk . . . '-0000000' 

where n.jk... is the number of first level cells in the ith 

second level cell, in the kth third level cell, etc. 

etc. to 

0+296363 *n. . . s '-0000000' 

0+2 96363 ' 1 ' - 0000000 * 

0+3 O6363 ' data ' -0000000 ' » 

As a check: £11 jk... = Number of observations 
Sn.jk... = Number of nijk. ..'s 
&i..k... = Number of n. ik«..'s 
b. Example of Input: 

The data from table 10.1U.1 of G. Snedecor. Statisti cal 
Methods. Iowa State College Press, 5th Ed. 1957. gives 
Calcium concentrations in leaves of turnip plants in units 
of % The data are: 



Plant 1 2 

Leaf 12312 



3 k 

23123 



Deter. 1 3.28 3.52 2.88 2.U6 I.87 2. 19 2.77 3.7U 2 55 ^78 km * xi 
2 3.09 3-W 2.80 2.1* 1.92 2.19 2.66 I'M I'M 3!8? h'2 Ifi 

The data tape would contain: 

Calcium Concentration in Leaves of Turnip Plants (<£,)' 1000 t-inv^ </< n? 

0+295000'2'2'2'2»2'2'2'2'2'2»2'2'-0000000« C^ 1000 times /o CR 

0+296363 '3 '3 '3 '3 '-ooooooo' cr " 

O+296363 'If'-OOOOOOO' CR 
0+296363* 1* -0000000' CR 

0+306363 , 3280'3090'3520'3480'. ..1io70'^o«3310'3310»-ooooooo' • 

?S3 ^ ft T° ? *?" min t tlons fo ' ea <* of the twelve leaves the 

of the Sur »wf ^^ S \ ,**' th6re 8Xe three leaves in ea <* 
of the four plants, hence iota- 3's are entered. There are four plants 

^JS, e SS ae f ' ^ ence l he k is entered next - ae in P ut of ^ 

numbers within levels must always terminate with 0+2 96363' 1'- 0000000'. 
Sth^f ?^^ d f^ ia1 ^ ^ginning with the determinations 

2JS? 2J* ^ ^ af °f f e firSt Plant Proceeding on to the determin- 
ations within the second le £ .f of the first plant, etc. 
d. Ouput: ' 

lm SntSns 1 f rd ° UPUt (n0t ***** thG ChangeS ° f 3 * h) for this exara P le 

a. Means: level 0. These are the mean determinations for each 

I ' I I Uld be noted that there BXe k gro^Ps of 3 means 
each. Each mean within a group is the mean for a leaf, and 
each group corresponds to a plant. 

b. Means: level 1. These four means are for the three leaves with- 
in the four plants and therefore they are the plant mean 
determinations . 

c. Means: Level 2. This mean is for the four plants in the 
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experiment (the grand mean), 
d. Analysis of Variance 

The Level Mean Square is for Determinations ^ri.thin Leaves. 
The Level 1 Mean Square is for Leaves vithin Plants. 
The Level 2 Mean Square is for Among Plants. 
2. The alternate output gives the mean square (variance) of determin- 
ations vithin each leaf printed after the leaf mean (Means: Level 0) 
in the alternate color. 
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